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Abstract. 

The experimentally observed amazing dependence of a critical magnetic field B c of a full 
field-induced spin polarization and a spin susceptibility x* °f a two-dimensional electron 
gas on the electron density can be explained by screening of the Coulomb potential. The 
possibility of spontaneous full spin polarization expected at lower electron densities also 
crucially depends on screening. 



1 Introduction 

Among spin related phenomena one of the most intriguing is a spontaneous spin 
polarization (magnetization) of a two-dimensional electron gas (2DEG). 

In absence of direct methods to measure a 2DEG magnetization only indirect 
evidence is so far available. The experiments of Shashkin et al. [1] revealed almost 
a linear dependence of a critical magnetic field B c of full spin polarization versus 
2DEG density n in Si inversion layers. The full polarization state was associated 
with the onset of a saturation of a magnetoresistance in a parallel magnetic field. 
The extrapolation of the experimental curve to B c =0 gave a non-zero value of n. It 
looked like a manifestation of a possibility of spontaneous spin polarization at lower 
electron density. However, the recent breakthrough to a more dilute 2DEG allowed 
to Pudalov et al. [2] to discover a substantial declination of B c (n) curve from the 
conjecture of Ref. [1]. A possibility of spontaneous spin polarization remains thus 
still unclear [3,4]. 

The main goal of present communication is to elucidate the impact of screening 
on a 2DEG spin polarization. Just screening results in the behavior of B c (n) and 
similar to that in the experiment. 

2 Model of induced and spontaneous spin polarization 

We adopt a simplified model for the two-dimensional Fourier transform of the 
screened Coulomb potential 
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where k is the permittivity, e is the elementary charge, and A is the model screening 
length. The latter could be defined by the 2DEG itself or by adjacent electrodes. 

The polarization degree r\ = \n+ — n_|/n where n+ and n_ are spin-up and 
spin-down electron densities, respectively, while n = (n + + n_) being the total 
electron concentration is derived via minimization of a total energy Etot including 
kinetic, interaction (exchange, correlation), and Zccman contributions. Here the 
direct Coulomb interaction is omitted as it does not depend on a spin polarization. 
We employ a quite common model 

ttH 2 1 
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where / is an interaction term, nib is the mean electron band mass, g v is the valley 
degeneracy, go is the free electron g- factor, /is is the Bohr magneton, and B is 
a magnetic field. The valley degeneracy g v equals 4 for the 2DEG in a silicon 
structur. 

For the non-interacting electron gas {1=0) the minimization of E to t with respect 
to a polarization degree 77 gives rise to a linear dependence of B c on n. Obviously, 
B c =0 when n—0. As for the spin susceptibility x* °^ v/E, it is constant and 
proportional to the product m&go=2TO&. For interacting electron gas it is convenient 
to express the spin susceptibility as x* oc m*g* where m*g* is a virtual product of 
electron mass and g-factor. 

For fairly small electron densities the interaction term in Eq. (2) becomes es- 
sential. 

Two plausible models of interaction could be used. The first one is that of 
Kohn-Sham (KS), widely employed for the interacting 2DEG description. The 
second one invokes the Ising (I) model widely exploited for spin lattice description. 
According for the KS model the interaction is 

I K s = I(n+)+I(n-). (3) 

Here the exchange and correlation belong to the electrons with parallel spins. If 
the interaction is associated with exchange, the values of I(n±) are supplied by 
integrals 

'^-~- s iLMLM^- kA <4) 

where the integrations are performed over the circle interior Q± corresponding to 
the Fermi wave vector Uf± = yjnn± (see Fig.l). In the limit Uf± » 1/A 

I(n±) ~ -(e 2 /nk F ±)n±. (5) 

i.e., as that for unscreened potential. In the opposite limit when kp± » 1/A 

I(n±) ~ -(e 2 A/ K )n|. (6) 

In the intermediate range of Uf± the interaction terms I(n±) may be numer- 
ically calculated via Eq. (4) for the screened Coulomb potential of Eq. (1). Evi- 
dently, the critical field B = B c should be associated with r) = 1. According to the 
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Figure 1: Two-dimensional momentum k-space of electrons in 2DEG; Uf+ and Uf- 
denote the corresponding Fermi wave vectors for spin-up and spin-down polarized 
electrons, respectively. 



Ising model the interaction term looks like 



h = Iks + .9, 




d 2 h 



d 2 k 2 



eV(k 1 -k 2 ). 



(7) 
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Here the interaction among electrons with opposite spin orientation is added. 
3 Spin polarization and conductivity 

The interplay of spin polarization and conductivity is worthy of being touched here 
as the experiments mentioned dealt with measurements of conductivity to detect 
the spin polarization. 
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Actually, a probable reason of the resistance growth and saturation observed 
in the experiments might be a lowered density of states (DOS) at the Fermi level 
caused by spin polarization. It could be also a likely cause of the metal-insulator 
transition if a spontaneous spin polarization (full or partial) may occur at zeroth 
magnetic field. Theoretically the lowered DOS at the Fermi level in the spin po- 
larized state was derived for 1DEG [5] to explain the observable declination of a 
quantum wire conductance from the conductance quantum, in particular, the so 
called "0.7 step". The existance of the lowered DOS ("pseudo-gap") at the Fermi 
level caused by spin polarization was also confirmed for 2DEG by calculations in 
[6]- 

Worth noting for the dependence of DOS deepening on the spin polarization 
degree to be substantial the Ising-like model for the exchange interaction should 
be employed instead of the conventional Kohn-Sham model. In detail this problem 
is discussed in [7]. 

4 Results and discussion 

The calculated dependences of the critical magnetic field B c on the electron sheet 
density n for the silicon structure are presented in Fig. 2. We used the Ising model 
for the exchange interaction given by Eq. (7) . The parameters of the structure are 
as follows. The electron mass equals to=0.36too where too is the free electron mass, 
k = (ksi + nsi0 2 )/^— 7-8 is the mean permittivity at the Si/Si02 interface, and 
the valley degeneracy is <j , „=4. The curve (1) corresponds to the model screening 
length A=13nm. The filled dots are taken from the experiment of Ref. [5]. The 
dashed line reproduces the conjecture of Ref. [1]. 

The calculated curve (1) in Fig. 2 to the best advantage coincides with the exper- 
imental points. It predicts a spontaneous spin polarization for n=0.4 10 11 cm~ 2 if 
the model screening length A fitted to experimental curve in Ref. [4] is actually as 
large as 13nm. However, the achieved electron densities are above 0.7 10 11 cm~ 2 
yet. For fairly small screening length A (for instance, if it approaches the Bohr 
radius) the screening results in both B c and n simultaneously come to zero and 
spontaneous spin polarization does not exist. This is illustrated by the curve (2) 
calculated for (A=5 nm). 

The calculated spin susceptibility \* °^ g*m* /2mi, against the parameter r s is 
plotted in Fig. 3. The ratio r s <~ (n) -1 / 2 /as characterizes the strength of Coulomb 
interaction with respect to the kinetic energy. The calculated curve in Fig. 3 also 
agrees well with experimental points (filled circles) because it is merely another 
presentation of results in Fig. 2. 

It should be emphasized that for unscreened Coulomb potential both models for 
the exchange interaction (KS and I) presented here inevitably lead to the existence 
of spontaneous spin polarization at sufficiently low electron density, that is, at high 
value of the parameter r s . 

The KS model was also used to fit to the experimental curve B c {n). The model 
screening length A turned out to be very large (about lOOnm) compared with 
that obtained for the Ising model. It could be said that the Ising model is more 
"resistive" against screening. 

As for the partial spontaneous spin polarization, i.e., r\ < 1, _B=0, in the frame 
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of the present models it could be hardly possible. Unlike to the induced spin 
polarization the spontaneous one has only an abrupt transition to a fully polarized 
state at sufficiently low electron density. 

In the experiments [1-4] the critical field B c , when the full spin polarization 
arose, was indirectly detected by the onset of a saturation of a magnetoresistancc. 
Therefore, it could be a discrepancy between the measured B c and the actual one. 
In our opinion, the critical field B c was underestimated in the experiments. 

Hopefully, the methods of a direct measurement of 2DEG magnetization will bc 
involved in the investigation. For example, in the recent experiments of Bagraev 
et al. [8] a nucleus magnetic resonance (NMR) of 29 Si atoms in the vicinity of a 
quantum wire was employed to detect its spin polarization. It was the first direct 
confirmation of the assumption that the so called "0.7 step" in a quantum wire 
conductance originates just in a spontaneous spin polarization. 



5 Conclusions 



The experimentally observed dependences of a critical magnetic field B c of a full 
field-induced spin polarization and a spin susceptibility \* °f a two-dimensional 
electron gas on the electron density are explained by screening of the Coulomb 
potential. The feasibility of spontaneous full spin polarization expected at lower 
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Figure 2: The calculated dependences /j,b-B c vs electron sheet density n (solid lines) 
for screening length A = I3nm (1) and A = 5nm (2). The dots correspond to the 
experiment of Ref. [5]. The dashed line reproduces the conjecture of Ref. [1]. 
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Figure 3: The calculated dependence of the parameter g*m*/2mb proportional to 
the spin susceptibility x* against the parameter r s for screening length A=13 nm. 
The dots correspond to the experimental data from Ref. [2] 

electron densities also crucially depends on this screening. 
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